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ABSTRACT 



< 



We present results from a near-infrared A'-band imaging survey of 59 elliptical 
galaxies in five nearby clusters. We measure photometric parameters for each galaxy 
using surface photometry and draw velocity dispersions from the literature. Three 
(3J[) observables define a near-infrared Fundamental Plane (FP) of elliptical galaxies with 

R e cx (T 1 - 44 ± - 04 x;~ - 79±0 - 04 . The scatter in the near-infrared relation is small at 16.5% 
in distance, which is equivalent to, or less than, the scatter of the optical FP. We 
£^ suggest that the small deviation of the near-infrared FP relation from the optical FP 

is due to the reduction of metallicity effects in the near-infrared bandpass. While the 
f — small scatter of the optical FP could be consistent with compensating effects of age and 

metallicity, the similarly small scatter of the near-infrared FP is nearly independent of 
00 metallicity and hence places a strong constraint on possible age spreads among elliptical 

galaxies at every point along the FP. We suggest that the departure of the near-infrared 
Q\ FP from the pure virial form R e cx cr 2 !!" 1 , and the corresponding observed relation 

(M I L) cx ji^f°- 16 ±°- 01 5 ma y be explained by slight systematic departures of the structure 
Qh and dynamics of elliptical galaxies from a homology. 
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1. Introduction 

Elliptical galaxies form a two-parameter family, a two-dimensional Fundamental Plane (FP) 
within the multi- dimensional parameter space of several observables (Djorgovski & Davis 1987; 
Dressier et al. 1987; see also Kormendy & Djorgovski 1989, Djorgovski 1992, and references 
therein). Typically, one considers a parameter space of the half-light radius r e , the mean surface 
brightness enclosed within that radius (/i) e , and the central velocity dispersion a. The FP 
correlations have a small scatter (of order 0.1 dex in logr e ), which makes them potentially useful 
both as distance indicators and as a strong constraint on galaxy formation models. Explanation 
has been offered as to the origin of the FP, in particular its apparent deviation from the pure 
virial relation, primarily by invoking stellar populations variations among galaxies along the FP 
(Djorgovski 1988; Djorgovski & Santiago 1993, hereafter DS93; Renzini & Ciotti 1993). Substantial 
fine-tuning of variations in the IMF or mass-to-light ratio are required to remain consistent with 
the thinness at every location along the FP. Furthermore, the allowed range of mass-to-light ratio 
from one end of the FP to the other is small (Djorgovski & Davis 1987; Faber et al. 1987). 

Variations in met alii city along the FP could give rise to the tilt of the FP in optical 
wavelengths, as line-blanketing effects in the metal-rich giant ellipticals would make their (/i) e 
fainter. Observations of the FP in the near-infrared could shed light on the origin of the 
FP by probing a wavelength region with a different sensitivity to stellar populations effects. 
Observations in the li'-band window (A = 2.2/im) will trace old stellar populations and the 
bolometric luminosity, as this bandpass is far less met alii city- sensitive than optical wavelengths. 
Early results on the near-infrared FP (Recillas-Cruz et al. 1990, 1991; DS93) have all relied 
on photoelectric aperture photometry, whereas the optical FP comparison data utilized detailed 
surface photometry. The rapid development in near-infrared imaging detectors has only recently 
made a near-infrared imaging survey of this scale feasible. 

The observations described in this Letter are preliminary results for part of an ongoing 
program to observe more than 250 early-type galaxies in eleven clusters. The complete details 
of the observational methods, reductions, calibrations, and comparisons with photometry from 
the literature will be given in a subsequent paper (Pahre et al. 1995b). These results follow 
preliminary results presented for the Virgo (Pahre et al. 1995a) and Perseus (Pahre et al. 1994) 
clusters individually, and for the entire five cluster sample (Djorgovski, Pahre, & de Carvalho 
1995) which includes the Virgo, Perseus, Coma, Abell 194, and Abell 2634 clusters. 

2. Observations 

Nine galaxies in the Virgo cluster were observed in the _fif s -band (2.16/im) at the 2.5 m du Pont 
Telescope at the Las Campanas Observatory in 1993 March, and are described elsewhere (Pahre & 
Mould 1994). The observations of the other four clusters were taken in the _fif s -band between 1994 
October and 1995 March using a new near-infrared camera (Murphy et al. 1995) on the 1.52 m 
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telescope at Palomar Observatory. This instrument is based on a NICMOS-3, 256 X 256 pixel 
array reimaged at 1:1 at the f/8.75 cassegrain focus, producing a 0.620" projected pixel size and 
a 159" field- of- view. The galaxy lists were mostly chosen from those of Lucey & Carter (1988) 
and Lucey et al. (1991a, b), except for the Perseus cluster, which was selected to include all 19 
galaxies in Faber et al. (1989). Each galaxy was observed for 15 X 60 s on-source, separated into 
five different telescope pointings, and the sky was estimated from a similar total number of frames 
taken before and after each galaxy. All observations were taken during photometric conditions, 
with a mixture of the UKIRT faint standards (Casali & Hawarden 1992) and the new HST faint 
standards (S. E. Persson, private communication) being used for calibration. Usually between 
5 and 15 standards were observed on each photometric night, which have a typical rms scatter 
less than 0.02 mag. The typical seeing was < 1.5", with some observations having significantly 
better seeing than could be fully-sampled. We have compared aperture magnitudes measured for 
the Virgo and Coma galaxies with Persson, Frogel, & Aaronson (1979), and find mean offsets of 
K ur — Kpfa = 0.00 ± 0.01 mag (Virgo) and +0.02 ± 0.02 mag (Coma), with an rms scatter of 
0.05 and 0.04 mag, respectively. 

We have constructed surface brightness profiles for each galaxy by finding the best-fitting 
elliptical isophote fj, in integer pixel steps in the semimajor axis a using the ELLIPSE task in the 
STSDAS package of the IRAF image reduction software. Profiles at a > 10 pixels were rebinned 
to improve the SNR as described in Djorgovski (1985), and then corrections to the photometry for 
surface brightness (SB) dimming as (1 + z)~ 4 , galactic extinction (taking Ak = 0.085A# and the 
median Ab estimate from Faber et al. 1989 for each cluster), and the i-correction of +3.32; for 
the li'-band (Persson et al. 1979). The mean SB internal to each isophote (/i) was calculated by 
integrating the profile at each semimajor axis length. Half-light effective radii r e and the mean SB 
interior to that radius (/i) e were measured using fits of fj, against a 1 / 4 under the assumption that 
all profiles followed a de Vaucouleur form. The apparent half-light radii r e were then transformed 
to metric size R e by assuming Hq = 75 km s _1 . We have drawn velocity dispersion a data from 
the literature, giving equal weights to the work by Lucey & Carter (1988) or Lucey et al. (1991a, b) 
and a mean of all other measurements from the literature (as compiled by McElroy 1994, private 
communication). We note that using the Faber et al. (1989) velocity dispersions for the Coma 
cluster causes no significant change in the FP parameters described in §3. We have excluded all 
galaxies with a < 100 km s _1 or r e < 2", thereby avoiding issues of dwarf galaxies and seeing 
corrections. The galaxy sample is thus 59 galaxies out of the initial list of 83. 

3. Constructing the Near-Infrared Fundamental Plane 

Using the measurements of r e , (/i) e , and a described in §2, we have used bivariate least-squares 
analysis to find the optimal correlation between these FP variables. The methodology is as found 
in previous studies (Djorgovski & Davis 1987; de Carvalho & Djorgovski 1989; de Carvalho & 
Djorgovski 1992). In order to construct the FP itself, we first find the "mixing" value b in the 
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quantity logu + b(fj,) e which has the highest correlation with logr e , and estimate the errors in b 
from the "bootstrap" method. We then measure the slope A and intercept C for the FP relation 

log r e (arcsec) = Aloga - 0AB{n) e + C, (1) 

where B = —2.5Ab. The parameters A and C are measured using a linear least-squares relation 
by minimizing the variance orthogonal to the best-fitting line. We have measured the parameters 
A, B, and C for all five clusters separately, and find that the values of A and B are consistent 
within the uncertainties. Taking the median values of A and B, we then measured the median 
intercept C" for each cluster, and used this value of C" to transform each galaxy onto a metric 
scale (under the assumption that the Coma cluster has no peculiar velocity and Ho = 75 km s _1 ). 
We then measured the value for A and B for the entire sample simultaneously in metric units for 
i? e , and find that the entire dataset is best-fit by the li'-band FP relation with A = 1.435 ± 0.040, 
B = —0.79 ± 0.04. This result is consistent with the median result from the individual clusters, 
and corresponds to the scaling relation r e cx ( 7 1 - 44 ± - 04 x;~ - 79±0 - 04 (where S e is the mean surface 
brightness within r e expressed in linear units, and (/i) e is the same quantity expressed in 
mag arcsec -2 ). The results for the individual clusters are detailed in Table 1 and displayed in 
Figure 1. 

We have calculated the quartile-estimated scatter in log R e to be Qa = 0.072 dex for the entire 
dataset, which corresponds to an error in distance of 16.5% per galaxy. This dispersion about 
the near-infrared FP is equivalent to, or smaller than, the 0.088 dex (Djorgovski & Davis 1987), 
0.07 dex (de Carvalho & Djorgovski 1992, for the Faber et ai. 1989 cluster sample), or 0.075 dex 
(Lucey et ai. 1991b) found for the optical FP relation. Measurement errors are typically 0.04 dex 
in logr e , 0.017 mag in fitting (/i) e , 0.02 dex due to sky-subtraction uncertainties, and 0.03 dex in 
log it, which suggests that much of the observed 0.072 dex thickness of the near-infrared FP can 
be accounted for by the 0.063 dex observational errors, resulting in an upper limit of 0.035 dex for 
the intrinsic thickness of the near-infrared FP. 

Our results for the near-infrared FP are consistent with previous photoelectric photometry 
results. For the li'-band, DS93 found A = 1.5 ± 0.1 and B = —0.78 ± 0.17 for the r e measurements 
of Faber et ai. (1989) combined with the colors of Persson et ai. (1979). Recillas-Cruz et ai. (1990, 
1991), found A = 1.69 ± 0.11 and B = -0.74 ± 0.05 for the Coma cluster and A = 1.55 ± 0.19 and 
B = —0.81 ± 0.10 for the Virgo cluster using their own li'-band measurements. 

We also investigated the Kormendy relation (see Kormendy & Djorgovski 1989 and references 
therein) of logi? e (kpc) = k(/2) e + /, which is an oblique projection of the FP. We assume the 
relative distances derived from the FP relation as described above in converting r e to R e . We use 
a least-squares minimization of the orthogonal distance from a linear relation between these two 
variables. We fit each cluster independently and find that they are all well-fit by values for k that 
are consistent within one-sigma of the mean (and median) value for k, therefore we can fit the 
entire sample to determine the mean relation. The best fit for the five clusters is k = 0.332 ± 0.023 
and / = —5.090 ± 0.022, corresponding to the scaling relation R e oc j- - 83 ± - 06 _ This relation has a 
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quartile-estimated dispersion of Qa = 0.173, corresponding to an uncertainty of 40% in distance 
per galaxy. This best-fitting Kormendy relation is also plotted in Figure 1. The near-infrared 
relation can be compared to the optical r-band relation of R e oc I~ ' 85 (Kormendy & Djorgovski 
1989). The near-infrared FP thus shows a 60% reduction in scatter over the near-infrared 
Kormendy relation. 

4. Discussion and Conclusions 

The best-fit to Eq. (1) of the near-infrared FP relation can be compared with its optical 
counterpart. The F-band relation is A = 1.24 ± 0.08 and B = —0.9 ± 0.1 (Lucey, Bower, & 
Ellis 1991, taking the mean of the Coma and Virgo cluster results) or A = 1.25 ± 0.07 and 
B = —0.80 ± 0.03 (the cluster sample of Faber et ai. 1989, in the analysis of de Carvalho & 
Djorgovski 1992). The variation in the exponent of the FP (i.e. the R e oc a A term) with bandpass 
was described by DS93, with the slope of the FP becoming steeper (A increases) with increasing 
wavelength. We find that A has changed by +0.19 ± 0.06 dex from the F-band to the li'-band, 
consistent with the DS93 result of +0.29 ± 0.11. 

The relationship in Eq. (1) departs, however, from that expected for a pure virial theorem, 
if light directly traces mass. The virial form of the FP is R e cx cr 2 !!" 1 , so the near-infrared FP 
formally differs from the virial relation at high confidence for the A parameter. Faber et ai. 
(1987) and Djorgovski (1988) argue that if elliptical galaxies form a homologous family, then the 
mass-to-light ratio enclosed by R e is a function of luminosity according to the scaling relation 
(M I L) e cx M a , where a = (2 — A) / (2 + A). [We note that the alternate form found in some 
studies is (M/L) e cx L%, where [i = (2 - A)/2A and a = /3/(l + fi).] The optical studies listed 
above result in a = 0.23 ± 0.02 in the F-band. The near-infrared result is a = 0.16 ± 0.01, which 
is inconsistent with zero to high significance. Since the near-infrared sample described in this 
Letter spans ~ 4 mag in luminosity, the least luminous galaxies have a mass-to-light ratio only 
50% less than that of the largest galaxies. Our near-infrared result is also consistent with DS93 
and Recillas-Cruz et ai. (1990, 1991) who found a = 0.14 ± 0.03 and 0.13 ± 0.05, respectively. 
This dependence of mass-to-light ratio on mass is plotted in terms of observables in Figure 2, 
although we note (as does Djorgovski 1988) that this figure shows an oblique projection of the FP, 
prone to large cumulative and correlated errors. Nonetheless, the FP-derived value of a = 0.16 is 
consistent with the data of Figure 2. 

Explanations of the departure of the optical relation from zero often focused on line-blanketing 
effects (Faber et ai. 1987) as there is a mass-metallicity relation along the elliptical galaxy 
sequence. But while the slope A of the near-infrared FP could be explained by stellar populations, 
the small scatter of the relation would be more challenging. For example, by assuming a 
uniform mapping between the color-magnitude relation and the FP, we find that the (V — K) 
color-magnitude relation for ellipticals (Bower, Lucey, & Ellis 1992) predicts a change in (V — K) 
of ~ 0.32 mag for the 4 mag range in Kj of our observations. In the Worthey (1994) stellar 
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populations models for solar met alii city and a 12 Gyr age, the color-magnitude relation would 
require only a change of 0.02-0.04 dex in (M/L)boI between the two extremes of the FP, which 
is less than the 0.3 dex required in the li'-band, thus the tilt of the FP cannot be explained by 
metallicity effects alone. A variation in age between 8 Gyr and 17 Gyr along the sequence would 
produce a 0.25 dex variation in (M / L)b u but this would require a synchronicity in all clusters 
in the formation of elliptical galaxies of a given mass, since the thickness of the FP at any point 
is only < 5% of its span in age. Another explanation for the variation in (M/L) along the FP 
could be the velocity anisotropy, but DS93 have used correlations utilizing velocity anisotropy 
and showed that they follow similar scaling relations to those described above; they conclude 
that velocity anisotropy could not explain all of the (M/L) dependence on mass. IMF variations 
may be partly responsible (DS93), but they require a fine-tuning which appears unlikely (Renzini 
& Ciotti 1993); the same may be said for a systematic variation in the dark matter content of 
elliptical galaxies. 

Under the assumption that elliptical galaxies form a homology (in both their spatial and 
kinematic distributions), there is a unique mapping between the observed central velocity 
dispersion and the mean kinetic energy per unit mass, and between the de Vaucouleur radius and 
the mean harmonic radius (Djorgovski, de Carvalho, & Han 1988). It is well known, however, 
that more luminous ellipticals tend to be more anisotropic (Davies et al. 1983) and to have 
shallower surface brightness profiles (Schombert 1986; Caon, Capaccioli, & D'Onofrio 1993). In 
addition, in A-body simulations there appear to be velocity dispersion aperture effects which 
entail a nonhomologous velocity distribution at any fiducial radius for galaxies of the same mass 
(Capelato, de Carvalho, & Carlberg 1995). Simple analytical models of elliptical galaxies can be 
constructed (Hjorth & Madsen 1995), which would reproduce the necessary surface brightness 
profile departures from the homology (Burkert 1993). The observations in this Letter that a ^ in 
the li'-band, suggests that the origin of much of the tilt of the FP (relative to constant (M/L) and 
the pure virial theorem) can be explained by systematic deviations of the structure and dynamics 
of ellipticals from a homology. In other words, the mappings described above may be changing 
along the sequence of elliptical galaxies. It is thus possible that the slope of the near-infrared FP 
can be explained primarily by dynamical arguments. In such a scenario in which slight departures 
from a homology are present within the family of elliptical galaxies, the small thickness of the FP 
might place strong constraints on galaxy formation (Capelato et al. 1995). More clues about this 
could be afforded by simulations spanning a wider range of galaxy properties, or by observations 
of the FP at higher redshifts. 
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Table 1. Individual Cluster Solutions for the Near-Infrared Fundamental Plane 



Cluster 


N 


A 


± 


B 


± 


C 


± 


Qo 


C 


Qo 


Abell 194 


13 


1.12 


0.10 


-0.76 


0.10 


-6.759 


0.021 


0.070 


-7.679 


0.068 


Abell 2634 


10 


1.44 


0.06 


-0.75 


0.13 


-6.696 


0.018 


0.051 


-8.037 


0.055 


Perseus 


16 


1.67 


0.17 


-0.79 


0.12 


-8.234 


0.024 


0.085 


-7.694 


0.087 


Coma 


12 


1.29 


0.08 


-0.71 


0.06 


-6.993 


0.016 


0.053 


-7.814 


0.065 


Virgo 


8 


1.62 


0.06 


-0.93 


0.05 


-8.436 


0.012 


0.036 


-7.083 


0.055 
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Fig. 1. — (a) The near-infrared FP of elliptical galaxies and (b) the Kormendy relation (an oblique 
projection of the FP) for the sample of 59 galaxies. The scatter is clearly reduced by the introduction 
of the central velocity dispersion term. The total scatter from the FP relation in units of log R e is 
given as Qa = 0.072, corresponding to a 16.5% uncertainty in distance per galaxy. Median error 
bars are also shown in each panel. 
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Fig. 2. — The variation of mass-to-light ratio (ordinate plus a constant) with mass (abscissa plus a 
constant) for the 59 galaxies in the sample. This diagram by itself is of limited utility in constraining 
a, as the errors in both axes are both cumulative and correlated. The data are well-represented, 
however, by the FP-derived a = 0.16, with a residual scatter of Qa = 0.11 dex about this line. 



